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QUASAR OPTICAL VARIABILITY AND BLACK HOLE MASS 
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RESUMEN 

Favor de proporcionar un resumen en espanol. If you cannot provide a Spanish abstract, the 
editors will do this. In order to investigate the dependence of quasar optical-UV variability on fundamental 
physical parameters like black hole mass, we have matched quasars from the QUESTl variability survey with 
broad-lined objects from the SDSS. Black hole masses and bolometric luminosities are estimated from Sloan 
spectra, and variability amplitudes from the QUESTl light curves. Long-term variability amplitudes (rest- 
frame time scales 0.5-2 yrs) are found to correlate with black hole mass at the 99% significance level or better. 
This means that quasars with larger black hole masses have larger percentage flux variations. Partial rank 
correlation analysis shows that the correlation cannot explained by obvious selection effects inherent to flux- 
limited samples. Wc discuss whether the correlation is a manifestation of a relation between BH mass and 
accretion disk thermal time scales, or if it is due to changes in the optical depth of the accretion disk with 
black hole mass. Perhaps the most likely explanation is that the more massive black holes are starving, and 
produce larger flux variations because they do not have a steady inflow of gaseous fuel. 

ABSTRACT 

In order to investigate the dependence of quasar optical-UV variability on fundamental physical parameters 
like black hole mass, we have matched quasars from the QUESTl variability survey with broad-lined objects 
from the SDSS. Black hole masses and bolometric luminosities are estimated from Sloan spectra, and variability 
amplitudes from the QUESTl light curves. Long-term variability amplitudes (rest-frame time scales 0.5-2 yrs) 
are found to correlate with black hole mass at the 99% significance level or better. This means that quasars 
with larger black hole masses have larger percentage flux variations. Partial rank correlation analysis shows 
that the correlation cannot explained by obvious selection effects inherent to flux-limited samples. We discuss 
whether the correlation is a manifestation of a relation between BH mass and accretion disk thermal time 
scales, or if it is due to changes in the optical depth of the accretion disk with black hole mass. Perhaps the 
most likely explanation is that the more massive black holes are starving, and produce larger flux variations 
because they do not have a steady inflow of gaseous fuel. 
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1. INTRODUCTION 

Variability is a characteristic property of quasars. 
It is non-periodic in nature and shows erratic 
behaviour. Optical-UV variability (in unbeamed 
quasar continua) occurs on time scales of days to 
decades, with fluctuations typically a few tenths of 
a magnitude, increasing at longer time scales. 

Variability is recognized as an important diag- 
nostic of the physical processes responsible for black 
hole activity. It is also important as the only method 
by which the smallest physical scales in AGN can be 
resolved. The behaviour of quasar optical variability 
with luminosity and redshift is well understood (see 
Vanden Berk et al. 2004 and references therein), but 
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the real cause of the fluctuations remain opaque. In 
order to shed light on the origin of quasar variabil- 
ity, relationships between variability and AGN fun- 
damental parameters, such as black hole mass and 
Eddington ratio, may be investigated. 

Here we describe a study that correlates quasar 
long-term variability with BH mass. Details are pro- 
vided by Wold, Brothcrton & Shang (2007). 

2. SAMPLE SELECTION 

Wc form a sample of 104 quasars by match- 
ing broad emission line objects in the SDSS DR2 
(Abazajian et al. 2004) with sources in the 200k 
Light Curve Catalogue of the QUESTl variability 
survey (Rengstorf et al. 2004). Redshifts were con- 
strained to z < 0.75 to keep the redshifted H/3 line 
within the spectral coverage of the SDSS, and to 
limit time dilation effects. 
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Fig. 1. BH mass and bolometric luminosity as a function 
of redshift for the SDSS-QUESTl sample. 

From single epoch Sloan spectra we estimate 
virial BH masses using the H/3 line width and the 
rest-frame quasar luminosity at 5100 A following 
the scaling relationships of Vestergaard & Peterson 
(2006). The bolometric luminosity of each quasar is 
evaluated as Lboi = 9 x AL^sioo a-nd the Eddington 
luminosity as iEdd ~ 1-51 x 10"^* Mbh/Mq ergs~^. 
Black hole mass and Lboi are plotted as a function 
of redshift in Figure [TJ 

2.1. Variability measurements 

Quasar variability is characterized by comput- 
ing the distribution of all possible variability am- 
plitudes (i.e. magnitude differences) on the quasar 



light curves, Ami 



rrij, where i < j. The 



variability of one quasar is described by the stan- 
dard deviation, mean, median and/or maximum of 
its distribution of Am (e.g. Givcon ct al. 1999). In 
Figure [5] we show the maximum variability ampli- 
tude for every quasar as a function of BH mass, and 
an increasing trend of variability with BH mass is 
seen. Sources displaying the largest variability am- 
plitudes have on average higher BH masses. The 
correlation is also illustrated in Figure [31 where it 
can be seen that most of the variability is detected 
at longer time scales. 

3. VARIABILITY-BH MASS CORRELATION 

Pearson correlation analysis shows that the corre- 
lation between BH mass and variability is significant 
at the w 3(7 level. The strongest correlation is mea- 
sured between BH mass and maximum variability 
(correlation coefficient 0.285 and a two-sided proba- 
bility of arising by chance of 0.3 %). Other measures 
of variability have qualitatively similar, but some- 
what less significant correlations with BH mass. 
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Fig. 2. Variability amplitude as a function of BH mass. 
The mean amplitude in five different BH mass bins is 
over-plotted as squares. 
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Fig. 3. BH mass versus the structure function, S(t), 
which measures how variability is distributed over rest- 
frame time scale r, here for time scales 1-100 and > 100 
days (up to «2 years for the SDSS-QUESTl sample). 



There are no significant correlations detected be- 
tween variability and Eddington ratio, iboi/^Edd, 
possibly because of the small sample size. The Ed- 
dington ratio is clearly an interesting parameter to 
correlate with variability as it measures whether 
variability is related to how effectively the quasar 
converts accreted matter into luminosity (assuming 
a standard optically thick accretion disk as described 
by e.g. Shakura & Sunyaev 1973). 
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3.1. Selection effects 

The variability-BH mass correlation could be 
caused by more primary correlations with redshift 
and luminosity as these are known to correlate with 
variability (e.g. Vanden Berk et al. 2004). For in- 
stance, higher-z quasars in the SDSS-QUESTl sam- 
ple are biased toward higher BH masses (Figure [T]). 
A correlation between variability and BH mass could 
therefore be produced if underlying correlations ex- 
isted between variability and z (or between variabil- 
ity and Lboi). But Spearman's partial correlation 
analysis shows that this is not very likely. E.g. the 
correlation coefficient between variability and BH 
mass, at constant z, is 0.203 with a probability of 
0.0035 of being caused by underlying trends with z. 
Similarly, Spearman's partial correlation coefficient 
between variability and BH mass, at constant Lboi 
is 0.254 with a probability of 2 x 10"'* of occurring 
because of underlying correlations with Lboi- 

Can contamination by residual host galaxy light 
in the spectra artificially enhance a correlation be- 
tween variability and BH mass? Residual host 
galaxy light may dilute the variability and would 
be a potential problem only for the lower-luminosity 
quasars. However, only very few spectra dis- 
play clear host galaxy features, and these were 
treated extra carefully during the spectral fitting. 
The variability-BH mass correlation is also present 
within a sub-sample of higher-luminosity quasars 
which are not expected to suffer from host contam- 
ination. Hence we conclude that although host con- 
tamination could artificially enhance the correlation, 
it is a weak effect and unlikely to be important. 

Finally, could the correlation be caused by time 
sampling effects? The best-sampled rest-frame time 
scale in the QUESTl survey is 1-2 years and this 
is also where most of the variability is detected 
(Rengstorf et al. 2006; Wold et al. 2007). Figure [3] 
shows that the variability at shorter time scales (1- 
100 days) corresponds to the noise level in the data 
(« 0.06 mag), but that the variability at longer time 
scales is intrinsic to the quasars. Hence there is 
a bias toward detecting intrinsic quasar variability 
only at longer time scales. If a relation between 
BH mass and variability time scale exists, this may 
lead to a preference for detecting variability only for 
quasars with higher BH masses. And quasars with 
lower-mass BHs would have fluctuations at or be- 
low the noise level, thus producing a variability-BH 
mass correlation. There are certainly arguments that 
characteristic variability time scale (such as accre- 
tion disk thermal time scale) may depend on BH 
mass (e.g. Collier & Peterson 2001), but this is still 



controversial and far from established. 

4. DISCUSSION 

We conclude that there is evidence for a correla- 
tion between the BH mass of a quasar and its optical 
long-term (0.5-2 yrs) variability properties. Quasars 
with more massive BHs tend to show larger variabil- 
ity amplitudes, i.e. quasars with larger BH masses 
have larger percentage flux variations. Our favoured 
explanation is that the correlation is real and not 
caused by selection effects. 

But what is its origin? Higher-mass BHs are 
more likely to be starving than those of lower mass. 
This is because higher-mass BHs must have swal- 
lowed more of their surroundings in order to have 
grown to their larger masses. Hence the larger flux 
variations for the more massive BHs could indicate 
that the BH is running out of gaseous fuel. Because 
there is not a steady inflow of gaseous fuel from the 
surroundings, larger fluctuations in the continuum is 
seen. 

The correlation could also be explained in terms 
of optical depth of the accretion disk. This can come 
about if the more massive BHs have accretion disks 
with larger optical depths than those of lower mass. 
A large optical depth indicates radiatively more ef- 
ficient accretion, hence any changes in the energy 
generation mechanism will be more easily transmit- 
ted through an optically thick disk. 

No significant correlations or anti-correlations 
could be detected between Eddington ratio and vari- 
ability, probably because redshift and luminosity ef- 
fects are difficult to disentangle in this sample. 

A larger sample, or a sample spanning a longer 
time base line, is needed to confirm the general ro- 
bustness of the variability-BH mass correlation. 
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